Two-channel Kondo physics in tunnel-coupled double quantum dots 
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We investigate theoretically the possibility of observing two-channel Kondo (2CK) physics in 
tunnel-coupled double quantum dots (TCDQDs), at both zero and finite magnetic fields; taking 
the two-impurity Anderson model (2AIM) as the basic TCDQD model, together with effective low- 
energy models arising from it by Schrieffer- Wolff transformations to second and third order in the 
tunnel couplings. The models are studied primarily using Wilson's numerical renormalization group. 
At zero-field our basic conclusion is that while 2CK physics arises in principle provided the system 
is sufficiently strongly-correlated, the temperature window over which it could be observed is much 
lower than is experimentally feasible. This finding disagrees with recent work on the problem, and 
we explain why. At finite field, we show that the quantum phase transition known to arise at zero- 
field in the two-impurity Kondo model (2IKM), with an essentially 2CK quantum critical point, 
persists at finite fields. This raises the prospect of access to 2CK physics by tuning a magnetic field, 
although preliminary investigation suggests this to be even less feasible than at zero field. 



I. INTRODUCTION 

Over the last decade or so, quantum dot devices 1 have 
become increasingly important testbeds for the realiza- 
tion and controlled experimental study of correlated- 
electron phenomena. The spin-i Kondo effect 2,3 is the 
classic example, in which at low temperatures the spin 
degree of freedom of the dot is screened as a result of tun- 
nel coupling to metallic leads. The rich physical behav- 
ior arising, in particular the strong many-body enhance- 
ment of the zero-bias conductance, 4-6 has stimulated the 
search for related phenomena in more complex device 
geometries: extensive work, both experimental and the- 
oretical, has uncovered a wide range of examples, includ- 
ing orbital and SU(4) Kondo effects, 7-11 underscreened 
Kondo behavior, 12-15 and several Kondo effects induced 
by an applied magnetic field, 16-20 to name but a few. 

In this paper we consider a tunnel-coupled double 
quantum dot (TCDQD). 21 The system consists of two 
locally correlated and mutually tunnel-coupled quantum 
dots, positioned in series between two metallic leads; 
and tunnel-coupled to them, such that current can flow 
through the system under a voltage bias applied to the 
leads. Experimentally, recent advances in nanofabrica- 
tion have enabled construction of such systems in both 
carbon nanotube 22-26 and semiconductor devices. 21 ' 27,28 

From a theorist's perspective, the canonical model de- 
scribing TCDQDs is the well known two-impurity Ander- 
son model (2AIM). 29-32 In a gate- voltage regime where 
each dot is effectively singly-occupied, the low-energy 
physics of the 2AIM is in turn embodied - to leading 
order in the tunnel-couplings under a Schrieffer- Wolff 
transformation 33 - in the two-impurity Kondo model 
(2IKM). 34-40 The physics of the 2IKM is immensely 
rich. 34-40 In particular, in the absence of an applied 
magnetic field, it is well known to contain a quantum 
phase transition, for which the quantum critical point 
is in essence a two-channel Kondo (2CK) fixed point 
(FP). 35, 36,39-41 That in turn raises the prospect of ob- 



serving 2CK physics in TCDQD systems, which recent 
theoretical work 42 has suggested to be potentially viable. 
This is a central issue considered in the present paper. 

The simplest exemplar of 2CK physics is the 2CK 
model, 43 consisting of a single spin-i coupled via an- 
tifcrromagnetic Kondo exchange to two metallic leads, 
which compete to Kondo-screen the spin and result in 
overscreening of it. 43 In consequence, the 2CK ground 
state is a non-Fermi liquid, characterized by the stable 
infrared 2CK FP and exhibiting exotic physical proper- 
ties such as a residual entropy of \ In 2 {k B = 1). 44 > 45 
The 2CK FP is however notoriously susceptible to desta- 
bilizing perturbations: 39 ' 46,47 inter-lead charge transfer in 
particular - as will always occur to some degree in a real 
device (and is inherently contained in the 2AIM) - is 
well known to destabilize the 2CK FP, 39 ' 46-50 rendering 
it unstable on the lowest temperature (T) scales. The 
system instead flows ultimately to a stable strong cou- 
pling (SC), Fermi liquid- like FP below some characteris- 
tic low-temperature Fermi liquid scale. 

For this reason, two-channel Kondo has experimentally 
been the most elusive of the various Kondo effects (we 
know of only one example 51 where it is believed to have 
been observed cleanly). Potential observation of it relies 
on the fact that if inter-lead charge transfer is sufficiently 
small, then a T-window can at least in principle exist over 
which the system flows close to the now-unstable 2CK FP 

- such that non-Fermi liquid behavior could be observed 

- before ultimately crossing over to the stable SC FP. 
The obvious questions then are: under what conditions 
does this arise, and are the resulting temperatures exper- 
imentally credible? 

These questions are considered in the present work 
where, for vanishing magnetic field in the first instance 
(sec. Ill), we study directly the 2AIM, together with 
the effective low-energy models derived from it (under 
Schrieffer- Wolff) to 2 nd and 3 rd order in the tunnel cou- 
plings; respectively the 2IKM, and a spin-model contain- 
ing the key effects of cotunneling inter-lead charge trans- 
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FIG. 1: Schematic of the 2AIM, as discussed in text. 



fer. We also consider for comparison the model studied 
in ref.42, in which solely direct inter-lead charge transfer 
is added to the 2IKM. The models themselves are speci- 
fied in sec. II, and the numerical renormalization group 
(NRG) method 52 " 54 (see ref. 55 for a review) is employed 
to study them, backed up by physical arguments. 

In sec. IV we consider these models when a non-zero 
magnetic field is applied to the dots. For the channel- 
symmetric 2IKM in particular, we show that its zero-field 
quantum phase transition is the terminal cndpoint of a 
line of transitions characterized by a 2CK FP. This raises 
the possibility that, in the presence of sufficiently small 
inter-lead charge transfer, 2CK physics might be acces- 
sible by tuning a magnetic field; which question is then 
considered. The paper ends with concluding remarks. 



II. MODELS 

We begin by specifying the models considered for 
TCDQDs, starting with the two-impurity Anderson 
model (2 AIM) 29-32 as the canonical model for such. 

A. Two-impurity Anderson model 

The 2 AIM is illustrated schematically in Fig. 1. It con- 
sists of two single-level dots (labeled v = 1,2, with level 
energies e v and on- level Coulomb replusion U) , mutually 
tunnel-coupled by a hopping matrix element t; and with 
each dot tunnel-coupled to a separate non-interacting 
metallic lead. The model Hamiltonian is 

#2 AIM = # leads + #dots + -Hhyb> (1) 

where for the double quantum dot itself, 
-Hdots =H e + Hjj + H t 

= ^2 ^n va + u X! ""t*w + 1 ( d i<r d 2a + H - c -) 

(2) 

with n va = S va d v(T the a ='f / ^-spin number operator 
for dot v. For the two equivalent leads (likewise denoted 



f = 1,2), 

Pleads = e k4ka C ^ka> ( 3 ) 

and we consider the standard case 2 of a flat-band lead 
with a (uniform) density of states per orbital of p = 
1/(2D), with D the half-bandwidth; denoting the total 
density of states by px = Np, with N (— > oo) the number 
of orbitals in a lead. The hybridization term coupling the 
dots and leads is 

Hhyb = J2 V » + H.C.) 

v^o (4) 

such that dot-^ is tunnel-coupled to lcad-i^ with matrix 
element V v \ and where 

fva — —fj^ c Ik<r (^) 

is the creation operator for the 'O'-orbital of the Wilson 
chain 2 ' 52 for lead v. Tunnel-coupling to lead-^ is embod- 
ied in the hybridization strength T v = ttptV^. Unless 
explicitly stated otherwise, we consider the case of sym- 
metric tunnel-coupling, Ti = T 2 = T; and the zero-bias 
Fermi level of the leads, Ep, is taken as the zero of energy. 

The effect of a magnetic field applied to the dots, which 
is considered in sec. IV, is encompassed by including 

Hb = -hS z (6) 

with S z — \ 5Z„(TVf — hvi) and h — gpsB; where B is 
the applied magnetic field and g the electron ^-factor. 

1. Symmetries 

If ei = e 2 and T x = T 2 = T (i.e. Vi = V 2 ), the 
model is 'left-right (LR) symmetric', meaning invariant 
under the transformation d\„ d 2a and cik CT <-> C2ko-- 
In addition, if ei = e 2 = —U/2 (but regardless of 
whether or not T\ = r 2 ), then the model is particle- 
hole (ph) symmetric, i.e. invariant under the ph trans- 
formation d\ a <-> d la ,d f 2a «• -d 2a , cJ k(T <-> Cl _ k(j and 

c 2k<T ^ — c 2-k<T- ^ n tms P a P er we consider explicitly 
the 2AIM at ph symmetry. 56 With LR symmetry also 
present, the full set of 'bare' parameters for the 2AIM is 
simply U/T, t/T and D/T. The bandwidth D is naturally 
taken to be the largest energy scale in the problem, and 
for our NRG calculations in practice we take D/T = 100. 

B. Schrieffer- Wolff transformations 

The 2AIM, allowing as it does for charge fluctuations 
on the dots, exhibits a rich range of behavior across its 
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full parameter space. Here we focus exclusively on the 
regime where each dot level is in essence singly occu- 
pied, as occurs for U 3> t, ri,r 2 in the ph symmetric 
systems considered. In this case a Schricffcr- Wolff (SW) 
transformation 33 may be used to obtain an effective low- 
energy model for the system. This involves perturbation 
theory in the tunnel couplings V\, V 2 and t, i.e. the per- 
turbing Hamiltonian is taken to be Hi — i?h y b + H t , and 
the only states of Hq = H — H\ retained are those in 
which the dots are singly occupied (with a local unity 
operator denoted 1'). 

This perturbation theory can in principle be carried 
out to any order in H\ . The leading non- vanishing contri- 
butions to the effective Hamiltonian arise to second order 
(specifically i'#i(e - Hq^PHiV, where P = 1 - V is 
a projector and eo is the ground state energy of Hq); and 
the effective low-energy model resulting from the second- 
order SW transformation on the 2AIM (neglecting re- 
tardation as usual 2 ) is the much-studied two- impurity 
Kondo model (2IKM). 34 - 38 ' 40 It consists of two spins- 
5, each coupled to a separate lead by antiferromagnetic 
(AF) Kondo couplings J\ and J 2 , and mutually cou- 
pled by an AF exchange coupling J - precluding as such 
charge transfer between the leads. The Hamiltonian is 

H2IKM — J\S\ ■ Sol + J2S2 • S 2 + -/Si • S 2 + -fflcads (7) 

where Sj, is a spin-i operator representing dot v = 1,2, 
and Squ is the spin-^ operator corresponding to the local 

spin density of lead v (S „ = fL* 7 **' f„a> with a 

the Pauli matrices). From the SW transformation, the 
parameters of the 2IKM are related to those of the 2AIM 
by pj v = 8T„/irU and J = U 2 /U. 

It is however obvious that a second-order SW trans- 
formation does not capture adequately the low-energy 
physics of the 2AIM, for it lacks the inter-lead charge- 
transfer processes that ensure the ground state of the 
2 AIM is always a Fermi liquid; and which are central in 
understanding the role of 2CK physics in TCDQDs. 

As mentioned above, SW to higher orders can be car- 
ried out, and to capture inter- lead cotunneling charge 
transfer one must go to third order. The additional 
third-order term arising from a SW transformation of 
the 2AIM, H 3 = i'H 1 (eo-Ho)- 1 PH 1 (e -H )- 1 PH 1 i', 
is given after lengthy calculation by 



-#3 = Vlr 



+ fLfia)Si • S 2 + 2A • (Si x S 2 ) 

(7 

(8) 



where V LR = {mV l V 2 )/U 2 = yJ{J.h,h)/U and A is a 
vector operator with components 
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So = S01 + So 2 (with a, /3, 7 € (x, y, z) and e a) 9 7 the Levi- 
Civita symbol). H 3 , in which charge transfer between 
the leads is mediated by the dot spins, is clearly a rather 
complicated object. It was obtained previously in Rcf. 
57, but subsequently neglected. In the following, we re- 
fer to the third-order effective low-energy model specified 
by -&21KM + H3 as the H3 model. 

It is important to emphasize that the charge-transfer 
processes in the 2AIM involve cotunneling, i.e. are me- 
diated by the dot spin degrees-of-freedom. In recent 
work 42,58 , the 2IKM with an additional direct lead- lead 
tunneling term was studied as a model for a TCDQD. 
The Hamiltonian considered was H 2 ikm + #3, with a 
charge transfer term 



H3 — Vlr y^(/lq/2g + flaf\a) 



(10) 



Note that A, which is self-adjoint and odd under 1 o 2 
exchange, is not a spin operator; its components satisfy- 
ing the commutation relations [A a , A@] = ze Q ,/3 7 5' Y , where 



where V' LR — jVlr- We refer to this as the H 3 model. 
It is not the correct effective low-energy model for the 
2AIM, and as such should not be expected to exhibit the 
same physics as the 2AIM even at low energies (indeed it 
does not ). We include it here purely for comparison to 
the 2AIM and H 3 models, to illustrate that adding the 
i?3 term to the 2IKM has a notably different effect to 
adding the H 3 term. 

We have now looked at all the models to be considered 
in this paper: the full 2AIM, the effective low-energy 
models for the 2AIM derived by SW transformation to 
2nd and 3rd order (the 2IKM and H 3 models respec- 
tively), and the H 3 model. In the following section we 
use the NIICJ 52 - 55 . 59 - 60 to obtain results for these models 
at zero field. We typically retain between 2000 and 4000 
states at each NRG iteration, and use an NRG discretiza- 
tion parameter A = 3. 



III. 2CK PHYSICS AT ZERO MAGNETIC 
FIELD 

We begin with a brief summary of the 2IKM at zero 
field, before considering the effect of inter-lead charge 
transfer as included in the 2AIM, H 3 and H 3 models. 



A. Two-impurity Kondo model 

It is well known 35-37,39 ' 40 that the LR-symmetric 
2IKM (Ji = J 2 ) exhibits a quantum phase transition 
(QPT) at a critical value J c of the inter-spin exchange 
J; with J c ~ 0(Tk), and 7k the Kondo scale of the sys- 
tem when the spins are decoupled, J = 0. The transition 
separates a local singlet (LS) phase arising for J > J c , in 
which the two spins bind to form a singlet, from a phase 
in which each spin is separately quenched by Kondo cou- 
pling to its attached lead (we refer to it as the 'Kondo 
singlet' (KS) phase). These phases are readily identified 
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from the phase shift, 5 e , in the even combination of con- 
duction channels/leads, which vanishes in the LS phase, 
and is ir/2 in the KS phase (see Refs.40,42,61 for details). 

The FP for the transition is distinct from those of 
the LS or KS phases, and corresponds to the 2CK 
FP, 35 ' 36,39 as known e.g. from conformal field theory 39-41 
and NRG 40 studies (albeit that the operator content 
and finite size spectrum of the critical FP differ slightly 
from the 2CK FP 40 ). On decreasing the temperature 
(T) /energy scale at J = J c , the system flows from a Lo- 
cal Moment (LM) FP - where the dot spins are effectively 
decoupled from each other and from the leads, with a cor- 
responding entropy Si mp = In 4 - to the critical 2CK FP 
characterized by Si mp = \ In 2, on the scale T <~ Tk (so 
that Tk is also in effect the two-channel Kondo scale). 

We comment in passing on the relation 35-38 ' 61 J c = 
q;Tk with a a constant, which we find from NRG calcu- 
lations indeed holds for sufficiently small pJ\ <C l- 62 The 
precise value of a naturally depends on how Tk (pertain- 
ing to J = 0) is defined; and in this there is freedom of 
choice. In practice we choose Tk to be the T for which 
^imp (Tk) = In 2, halfway between S'imp = hi 4 character- 
istic of the LM FP and the T = entropy S'imp = for the 
stable strong coupling FP in the KS phase. With this, the 
constant a rts 8. If instead we had chosen to define Tk as 
8/2.2 times the temperature for which S; mp = In 2, then 
a ~ 2.2; as often quoted in the literature. 35-38,61 There 
is however no fundamental distinction between these dif- 
ferent practical definitions of Tk- 

Finally, while the comments above refer to the LR- 
symmetric case, we add that the QPT is also known 63 ' 64 
to remain robust to J\ ^ J 2 , with a line of 2CK critical 
FPs in the (Ji, J2)-plane separating LS and KS phases; 
and a critical J c dependent on T^F* and the two 

distinct J = Kondo scales now arising. 



B. Effects of charge transfer 

We now look at the effect of adding inter-lead charge 
transfer processes to the 2IKM (focussing on the LR- 
symmetric case). These destroy the QPT occurring in 
the 2IKM, and with it the stability of the 2CK quantum 
critical point; the pristine transition being replaced by a 
continuous crossover between KS and LS ground states, 
characterized by a stable SC FP with Si mp (T = 0) = 0. 
This is well known for the 2AIM 48-50 and i^ 42 ' 47 models; 
and our NRG calculations indicate the same for the H3 
model (unsurprisingly, it being the effective low-energy 
model for the 2 AIM). 

Although the 2CK FP is rendered unstable by charge 
transfer, with decreasing T the system may first flow close 
to it on a scale T <~ Tk (as for the 2IKM) - evident e.g. 
in a characteristic ^ In 2 entropy plateau - before flowing 
to the stable SC FP on a low-energy Fermi liquid scale T* 
[which we calculate in practice from Si mp (T*) = jln2, 
halfway between the characteristic values for the 2CK FP 
and stable SC FP]. If this situation arises, then the 2CK 




FIG. 2: Simp vs pT (p = 1/(273)) from NRG for: 2AIM with 
U/T = 20 and t = U (solid line), and for the 2IKM (long- 
dashed), H3 (dotted) and H' 3 (short-dashed line) models with 
pJi — 0.093 and J — J c . The common Kondo scale Tk is 
indicated by an arrow. 



FP is effectively 'visible' at finite temperature, occurring 
over an appreciable T-window provided 

T* C T C Tk. (II) 

The obvious questions then are 42 : under what conditions 
does this behavior arise? And for the 2AIM in particular 
(as the canonical model for TCDQDs), docs it occur for 
experimentally realistic temperatures? 

To answer these questions we consider explicitly the 
T-dependence of the entropy Si mp (T). As for the 2IKM, 
Tk is the Kondo scale when the spins/dots are decou- 
pled, viz. J—Q for the spin-models T/3 and T/3, and t = 
for the 2AIM; with T K defined via Sim P (T K ) = In 2, as 
above (although in practice the resultant Tk differs in- 
significantly from that which can be read off e.g. fig. 2 
below from Si mp (TK) = In 2). To optimise the possibility 
of observing the 2CK FP at finite-T we follow ref. 42 and 
consider J = J c (~ Tk) for all spin models and t = t c for 
the 2AIM - where the models flow closest to the 2CK FP 
- chosen 42 in cither case so that the even channel phase 
shift S e = 7r/4. The phase shifts S e are themselves de- 
termined straightforwardly from the potential scattering 
on the even lead at the SC FP 54 (itself obtained by com- 
paring the NRG FP energy levels with those calculated 
separately from free even and odd conduction chains with 
equal and opposite potential scattering). 

Fig. 2 shows NRG results for S imp (T) for the 2AIM 
with U/T — 20 (solid line), with Tk as indicated. We 
wish to compare this to spin- models (H3, and 2IKM) 
with the same Tk, to which end we consider p,J\ = pj 2 = 
0.093 (which value differs somewhat from that given by 
the SW transformation, reflecting simply the fact that 
for U/T = 20 the SW asymptotics for pJi have not quite 
been reached). And for the spin models with charge 
transfer, we take Vlr = \/ J C J'[/U as in sec. II B. 

On decreasing T for the 2IKM with J = J c , the sys- 
tem naturally flows to the 2CK FP (Sj mp = \ In 2) for 
T ~ Tk, and remains there. For the T/3 model, with 
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(J-J c )/J c T/T k 



FIG. 3: For the same parameters as Fig. 2, the T = zero- 
bias conductance G c vs (J — Jc)/Jc for the 2AIM (solid line), 
the H3 model (long-dashed) and the model (short-dashed). 



FIG. 4: 5 imp (T) vs T/T K for the 2AIM (with t = t c ) for 
U/T = 20 (solid line), 30 (long dash), 40 (short dash), 50 
(dotted). A iln2 plateau appears as U/T is increased, indi- 
cating flow in the vicinity of the 2CK FP. 



direct inter-lead tunneling, there is a clear \ In 2 entropy 
plateau before flow to the stable SC FP (S[ mp — 0) on the 
scale T <~ T*; with T* < T K in this case, so that 2CK 
physics arises over an appreciable temperature window. 
The behavior of the 2AIM for U/T = 20 is quite different. 
After descent on the scale T ~ U from its trivial high-T 
limit S imp = In 16 to the LM FP with S* imp = In 4, the 
system flows directly to the SC FP on the scale T ~ Tk; 
with no hint of flow in the vicinity of the 2CK FP. The 
H 3 model is also seen to exhibit the same low-energy 
behavior. 

NRG may also be used to calculate the T = zero-bias 
conductance, G c . For LR-symmetric systems, 61 

G c = ^ sin2 (^ - S ) = 2 4- sm 2 (2S e ) (12) 
h h 

with <5 e ( ) the phase shift in the even (odd) channel, and 
the second equality follows at ph symmetry. Calcula- 
tion of 6 e thus gives G c directly. Fig. 3 shows G c vs 
[J - J c ]/J c for the 2AIM (where J = At 2 /U is taken), 
and the H 3 and H' 3 models; with J c <~ Tk as ever. For 
the H 3 model the halfwidth of this conductance peak (vs 
[J — J c ]/ J c ) is known 42 to be proportional to ^/T*/Tk; 
the fact that it is evidently <C 1 indicating the clear scale 
separation T* <C Tk seen from the entropies of fig. 2. G c 
for the 2AIM and H 3 models are similar, as expected (the 
differences again reflect that SW is asymptotically exact 
only as U/T — > oo). In these cases, by contrast, the con- 
ductance halfwidth is clearly 0(1). This too is consonant 
with the entropies shown in fig. 2, where for U/T = 20 
there is no flow in the vicinity of the 2CK FP, and as 
such Tk is the sole low-energy scale in the problem. 

We also point out here that our conductance results 
are in agreement both with previous work on the 
2 AIM 50 , and recent work on the H 3 model 42 : that the 
two models give different results is a natural conse- 
quence of the fact that they are not simply related by 
SW transformation. 

The 2CK FP, as manifest in a ^ In 2 entropy plateau, 



does not then occur for U/T = 20 in the 2AIM (or in 
its effective low-energy model, H 3 ). To observe it in the 
2AIM necessitates a larger U /T in order to suppress inter- 
lead cotunneling charge transfer (i.e. to reduce Vlr in 
the effective low-energy model), although this will also 
reduce T K itself (since 2 T K oc exp(-nU/{8T))). To illus- 
trate this, Fig. 4 shows S imp (T) vs T/T K for the 2 AIM, 
for various values of U/T. On increasing the interaction 
a ^ In 2 entropy plateau appears, indicating the opening 
of a temperature window in which the system flows in 
the vicinity of the 2CK FP, with a clear scale separation 
T* < T K for sufficiently large U/T, in practice U/T > 40. 
In this regime our numerics are consistent with the form 

and although we have not performed exhaustive calcula- 
tions our results indicate F(x) <~ bx 2 with b <~ 0(1) a 
constant. The behavior T* oc T^ 2 - also arises 42 in the H 3 
model. 65 Here it is found 42 that 

T* - b'( P Vl R ) 2 T K (14) 

with b' <~ 10 2 approximately constant, V' LR — \Vlr and 
(as above) (pVlr) 2 — J c (pJi) 2 /U . Since J c ~ T K itself, 
eq.14 gives T* cx T 2 . 

As mentioned above, increasing U/T in the 2AIM in 
order to access a reasonable T-window for 2CK physics, 
naturally has the effect of reducing Tk, the minimum 
temperature needed to access the 2CK regime. From 
fig. 4, U/T > 40 is in practice required for a reasonable 
window to arise. But for U/T = 40, our NRG results 
give T K /U - 4 x 10~ 9 . Taking U = 2 meV, as is 
experimentally typical 23,26 ' 28,66 , this corresponds to 
T K - 10~ 7 K - far lower than can be reached in exper- 
iment (a minimum of around lOmK). Even reducing 
U/T to 30, where a 5 In 2 'plateau' is just about visible 
in fig. 4, yields Tk ~ 10~ 5 K, some three orders of 
magnitude lower than what is experimentally feasible. 
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In contrast to previous estimates, 42 our conclusion is 
that it is highly unlikely 2CK physics could realistically 
be observed in tunnel-coupled DQDs. 

While we have focussed above on LR-symmetric sys- 
tems, the physics of the models considered remains qual- 
itatively the same when LR symmetry is broken: inter- 
lead charge transfer, whether of cotunneling or 'direct' 
form, destroys the QPT exhibited by the 2IKM. If we 
consider breaking LR symmetry by decreasing just one 
of the dot-lead couplings, then inter-lead charge trans- 
fer is suppressed, leading to a decrease in the T* scale. 
This raises the question whether 2CK physics might more 
readily be observed in a strongly asymmetric TCDQD 
device. That is not however the case. As for the LR- 
asymmetric 2IKM, the (experimentally relevant) 2AIM, 
if it flows at all in the vicinity of the 2CK FP, does so 
on the scale of min 

( t k ) > T k ) )- Decreasing e.g. T 2 at 
fixed Ti will thus not only reduce the T* scale, but also 

(2) 

the scale T K on which the system flows to the 2CK 
FP. Introducing LR-asymmetry in this way therefore de- 
creases the temperature at which 2CK physics might be 
observed. 



IV. 2CK PHYSICS AT FINITE MAGNETIC 
FIELD 

Having discussed the possibility of observing 2CK 
physics in a TCDQD at zero magnetic field, we now ask 
whether it might be possible to do so at finite field. Again 
we begin with the 2IKM, showing first that in the LR- 
symmetric case the quantum phase transition known to 
arise at zero-field is the terminal point of a line of QPTs 
accessed by tuning the magnetic field. 



A. Two-impurity Kondo model 

At zero field, the trivial atomic-limit of the 2IKM (i.e. 
eq.(7) with J\ = J 2 = 0) has a singlet ground state, 
15*) = -^(| t4-) — I It)) (representing the two dot spins), 
with a degenerate triplet at energy J above the ground 
state. Application of a magnetic field to the dot spins 
(eq. 6) lowers the energy of one triplet component, \T\) = 
I tt) (f° r h > 0), with a triplet-singlet energy difference 
E Tl - E s = (J - h). At a 'critical' field h = h c = J the 
ground state will thus be doubly degenerate, constituting 
as such as pseudospin- i comprised of \S) and |Ti). 

The energies of the remaining \Tq) and |T_i) triplet 
components lie at least J above the ground state. On 
coupling the dot spins to the leads, only the pseudospin- 
\ need therefore be retained in the low-energy mani- 
fold of dot states, provided J >• Tk ~ J c (or J 3> 

maxf^'.Tf) for Ji ? J 2 ); since Tk (defined as usual 
for J = 0) is the sole low-energy scale in the problem 
when J = 0. In othcrwords, the \T ) and |T_i) triplet 
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FIG. 5: NRG-determined phase diagram for LR-symmetric 
2IKM as a function of h and J, for pJi = pJ-2 = 0.093. J c 
(~ Tk) is the critical J for zero field (pj c ~ 1.63 x 10~ 5 here). 



components may be neglected provided J » J c , and only 
the pseudospin- ^ need be retained. Since this pseudospin 
can be flipped by the Kondo exchange on coupling to 
the leads, and since it is coupled to two leads, we thus 
expect 19 that two-channel Kondo physics should arise in 
the 2IKM with J ^> J c , at a critical field h = h c <~ J. 

The preceding physical argument may be put on firmer 
footing by deriving an effective low-energy model for the 
2IKM with J » J c and h ~ h c , retaining only the states 
which form the two components of the pseudospin (with 
local unity operator V = \S)(S\ + |T\><T\|); i.e. from 
H c s = I'^ikmI' to leading order, with H21KM in eq.7. 
This yields 67 the following effective low-energy model 



H c s = H' 



off 



O: 



+ 



(15) 



(in addition to H\ C ads, taken as read). Here r denotes 
the pseudospin- i, with components f z = i(|T 1 )(T 1 | — 



\S)(S\), 
with 



|Ti)(5| and f~ = (f+)+; and O = O x + 6 2 



0\ —jJiSq! + 

62 ={J - h)r z 



4J2S02 



(16) 



The first two lines of iJ ff , denoted H' cff , comprise a two- 
channel Kondo model with both channel anisotropy (for 
Ji 7^ J2), which is well known to be a relevant perturba- 
tion to the 2CK model; 39,46 ' 47 and with spin anisotropy, 
known to be irrelevant to the 2CK FP. 46 ' 68 

Since any channel anisotropy destroys the 2CK FP, 
consider first the LR-symmetric case J\ = J 2- Now we 
must consider the effect of O = 0\ + O2, each term of 
which is a relevant perturbation, 68 separately rendering 
the 2CK FP unstable (as we have also confirmed directly 
via NRG on eq.15). However for any given J, J\ and J 2 , 
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FIG. 6: S imp (T) vs pT for 2IKM with pJi = 0.093 and 
pj = 5 x 10~ 4 (~ 31pJ c 3> pJc) at the critical field h — h c 
(ph c — 4.55 x 10~ 4 or h c /J — 0.91), upon decreasing J2 from 
Ji: 2p(Ji - J 2 ) = (point-dash line), 10~ 6 (dotted), 10~ 5 
short dash), 10 -4 (long dash) and 10~ 3 (solid); correspond- 
ing respectively to [Ji — J2]/ J c = 0, 0.03, 0.3, 3 and 30. 




FIG. 7: S irnp (T) vs pT for different models with h ~ ft c . 
Shown for the 2IKM (short dashed line), H' 3 model (dotted 
line), and if 3 model (long dash); in all cases, p.J\ — p.Ji = 
0.093, pJ = 5 x 10~ 4 (> pJ c ), and Vlr = y/JJuh/U. For 
the 2AIM (solid line), we consider U/T = 20 and t/T = 0.71 
(corresponding to pJ = p[4t 2 /U] — 5 x 10 -4 ). 



the magnetic field h is a free parameter; which can then 
be tuned to ensure a vanishing coefficient of the relevant 
primary field associated with O, rendering it ineffective 
and the 2CK FP in consequence stable. This is the criti- 
cal field, h = h c . From the physical arguments above we 
expect h c ~ J (although not identically J, as is obvious 
from eq.16). 

The above arguments imply that for the LR-symmetric 
2IKM (Ji = J 2 ) with J > J c , there should be a QPT at 
a critical h c , with a 2CK critical FP; and hence a line of 
2CK critical FPs in the (h, J)-plane. We have confirmed 
this with NRG calculations on the 2IKM. Fig. 5 shows 
an illustrative phase diagram for fixed pJ± = pJ 2 as a 
function of (h/J c ,J/J c ) (recall that J c is the critical J 
for zero-field, with J c <~ Tk); the critical line h c (J/J c ) 
separating an h > h c phase which is continuously con- 
nected to the zero-field KS state, from that connected to 
the zero-field LS phase. Although the arguments given 
above apply strictly to J ^> J c , the transition is seen 
to extend continuously down to J = J c (as a simple ap- 
peal to continuity would suggest). As expected from the 
physical arguments above, h c ~ J indeed arises for suf- 
ficiently large J. Indeed near-linear behavior is seen in 
practice to set in for J/J c > 2 or so, and for J/J c ^> 1 
we have confirmed the linear form h c = —a + bJ where 
the gradient 6 — >• 1 (and a > is a constant). 

For Ji 7^ J 2 by contrast, H c ff (eq.15) is chan- 
nel anisotropic. The 2CK FP is consequently 
unstable, 39 ' 46,47 as likewise follows from the arguments 
above for the 2IKM at h = h c (at least for J ^> J c ). The 
system instead flows to a stable SC FP with <Si mp (T = 
0) = 0; flowing for sufficiently small | J\ — J 2 | from the 
2CK critical FP (S imp = \ In 2) to the stable SC FP, on 
a scale T* known from the 2CK model 46,69 to vanish as 
T* ~ (Ji — J 2 ) 2 . The validity of this picture has been 
confirmed by NRG and is illustrated in fig. 6, showing 
the T-dependence of S irnp (T) for the 2IKM at the criti- 



cal h = h c for fixed J and J\, upon decreasing J 2 from 
J 2 = Ji where the 2CK critical FP is stable. The low- 
energy scale T* - which in practice may be identified from 
<Simp(I*) = j In 2 - is immediately evident on increasing 
Ji — J 2 from zero; and analysis of the numerics indeed 
confirms it to vanish as T» ~ (Ji — J 2 ) 2 . 

We also note that the temperature scale (call it T^) 
on which the \ In 2 entropy plateau in fig. 6 is reached, 
is visibly lower than its counterpart shown in fig. 2 
for the 2IKM at zero field (with the same pJi), which 
is Tk oc exp(— 1/pJi). 2,70 This can be understood from 
the low-energy model H' eS (as appropriate to fig. 6), it 
being sufficient to consider the channel-symmetric case 
Ji = J 2 . This is a spin-anisotropic 2CK model, of form 
H' cti = E v =i,2[J*^' + y ± {§+T-+§^f+)] with (see 
eq.15) exchange couplings J z = Ji/2 and Ji = Ji/y/% 
(each less than Ji). And the characteristic low-energy 
scale for the model, T^,, on which temperature scale the 
i In 2 entropy is approached, is readily shown from per- 
turbative scaling 2 ' 70 to be oc cxp(— ^ -^-), whence 
T K «T K . 

This means, in othcrwords, that much lower temper- 
atures are needed to observe 2CK physics at the criti- 
cal point of the 2IKM at finite field than at zero field - 
scarcely a viable prospect in the light of sec. Ill B. 



B. Two-impurity Anderson model 

We now turn to the 2AIM at finite field. If t > t c 
(i.e. J = 4t 2 /U 3> J c ), we expect similar behavior to the 
2IKM: a sufficiently large magnetic field <~ h c renders the 
singlet and lowest triplet atomic-limit states degenerate, 
suggesting the possibility of 2CK physics at low tem- 
perature. As was the case in sec. Ill B however, whether 
two-channel Kondo can actually be observed in the 2AIM 
at finite field depends on the effect of charge transfer in- 
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herent to the model, since this destabilises the 2CK FP. 

Following the approach of sec. IIIB, we compare the 
behavior of the (channel-isotropic) 2AIM with two effec- 
tive low-energy models: the model in a magnetic field 
h <~ h c , and the result of a third-order Schrieffer- Wolff 
transformation of the 2AIM. The latter is found to be 
equivalent to the H 3 model at h ~ h c upon retaining only 
the \S) and \T{) dot states (viz H^ s = 1'(H 2 ikm+H 3 )1'), 
and is given explicitly by 

H'Jb = H cS - V LR [£ *(flf 2a + flf la )(r z I) 

<J 

(17) 

where a = +/— <^=^ t / 4-> an d -^eff is given by eq. 15. 

Fig. 7 shows results for Si mp (T), in which the pa- 
rameters are chosen to correspond to the 2AIM with 
U/T = 20, as considered earlier in fig. 2. First, for com- 
parison, the short-dashed line reproduces the 2IKM re- 
sult from fig. 6 for J\ = J2, where the 2CK critical 
FP is stable and hence 5 imp (T = 0) = ±ln2. Next, 
the dotted line in fig. 7 shows the behavior of the cor- 
responding iJg model (direct inter-lead charge transfer) 
with Vlr = \J JJ1J2/U. Since destabilizes the 2CK 
FP there is ultimately a crossover to the SC FP with 
vanishing residual entropy, but the crossover scale here 
is sufficiently small compared to that a ^ In 2 entropy 
plateau remains visible. The temperature window over 
which the plateau exists can be optimised by tuning h 
very slightly away from the critical h c of the correspond- 
ing 2IKM; this has already been performed for the dotted 
line in fig.7 (where (h — h c )/h c ~ 2 x 1CT 4 ). 

Turning now to the 2AIM and H 3 models, we have 
repeated the process of searching for the widest \ In 2 
entropy plateau by varying h around h c . However our 
calculations show no such plateau for the parameters con- 
sidered: the entropy always crosses directly from In 2 to 
zero, as illustrated by the solid and dashed lines in fig. 7 
(2 AIM and H 3 models, respectively). This is not surpris- 
ing, since it mirrors the h = behavior shown in fig. 2: 
while the model for U/T — 20 shows a distinct 2CK 
entropy plateau, the 2 AIM and H 3 models do not. 

In sec. IIIB we explained that the 2CK FP can indeed 
be observed at h = in the 2AIM, if one considers the 
model at a larger U/T > 40. It is thus natural to ask if the 
same is true for 2CK at finite field. We have undertaken 
preliminary calculations for U /T up to 50 (in practice the 
largest for which the calculations are feasible) in an effort 
to answer this question. In each case we have examined 
the entropy around h ~ h c <~ J c for signs of a | In 2 
plateau but, as illustrated in fig. 8, do not find any sign 
of 2CK behavior. This suggests that if 2CK is to be found 
in the 2AIM at finite field, it will arise only for U/T in 
excess of 50, for which the corresponding T' K will surely 
be out of range of experimental grasp. 
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FIG. 8: S imp (T) vs pT for the 2AIM with h ~ h c , t/T = 1 > 
tc/r, and U/T = 20 (solid line), 30 (long dash), 40 (short 
dash) and 50 (dotted). As U/T increases, the In 2 entropy 
plateau extends to lower T, but for the values of U /T consid- 
ered no \ In 2 plateau is found. 



V. CONCLUDING REMARKS 

In this paper we have examined the possibility of 
observing two-channel Kondo physics in tunnel-coupled 
DQDs. In the two-impurity Kondo model limit, such 
physics clearly arises at zero magnetic field near the 
2CK critical point, but in real quantum dots the effects 
of inter-lead charge transfer, which destabilise the 2CK 
fixed point, must of course be considered. 

While direct inter-lead hopping generally destroys the 
2CK physics on energy scales small compared to the 
Kondo scale, 42 thus providing a seemingly large window 
over which 2CK behavior should be observable, we have 
argued that cotunneling charge transfer processes - pro- 
ceeding via the dot spins, and arising naturally within 
the 2AIM - significantly reduce the likelihood of realis- 
ing the 2CK physics experimentally. 

A finite magnetic field opens up the possibility of a 
field-induced 2CK effect. For channel (LR-) symmet- 
ric systems, we showed that the quantum phase tran- 
sition of the zero-field 2IKM is the terminal point of a 
line of QPTs at finite field, the effective low energy criti- 
cal model at large h being a spin-anisotropic 2CK model. 
Again, however, the charge transfer processes present in 
real TCDQDs turn the line of QPTs into a crossover; in 
this case their effect is even more destructive, and we find 
no evidence of field-induced 2CK physics in the 2AIM on 
experimentally-realisable energy scales. 

It has been proposed that longer, even-numbered quan- 
tum dot chains might be good candidate systems for ob- 
serving 2CK physics. 63 Increasing the number of dots be- 
tween the leads suppresses charge transfer, 63 but it also 
leads to a decrease 64 in Ik, so that while the longer dot 
chain systems are more likely to flow to the 2CK FP 
(7k 3> T*), they are likely to do so at lower tempera- 
tures than for the two-dot case considered here. More 
work is needed to determine whether longer dot chains 
indeed offer a more promising route to accessing 2CK 
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physics. 
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